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Abstract 

Using conditional measurement on a beam splitter, we study the transformation 
of the quantum state of the signal mode within the concept of two-port non- 
unitary transformation. Allowing for arbitrary quantum states of both the input 
reference mode and the output reference mode on which the measurement is per- 
formed, we show that the non-unitary transformation operator can be given in a 
closed form by an s-ordered operator product, where the value of s is entirely de- 
termined by the absolute value of the beam splitter reflectance (or transmittance). 
The formalism generalizes previously obtained results that can be recovered by 
simple specification of the non-unitary transformation operator. As an applica- 
tion, we consider the generation of Schrodinger-cat-like states. An extension to 
mixed states and imperfect detection is outlined. 



1 Introduction 



Entanglement is one of the most striking features of quantum mechanics. Roughly 
speaking, a quantum state of a system composed of subsystems is said to be en- 
tangled, if it cannot be decomposed into a product of states of the subsystems and 
the correlation is nonclassical (note that there is no generally accepted definition 
of the degree of entanglement []1|). Recently applications of entangled quantum 
states in quantum information processing have been extensively discussed f^. 
Entangled quantum states also offer novel possibilities of quantum state engi- 
neering using conditional measurement. One of two entangled quantum objects 
is prepared in a desired state owing to the state reduction associated with an 
appropriate measurement on the other object. The quantum state of travelling 
optical modes can be entangled, e.g., by mixing the modes at an appropriately 
chosen multiport. The simplest example is the superposition of two modes by 
a beam splitter. Combination of beam splitters with measuring instruments in 
certain output channels may therefore be regarded as a promising way for engi- 
neering quantum states of travelling optical fields P, H, ^ H, ^ |], ^ |10], 0]. 

The action of a beam splitter as a lossless four-port device is commonly de- 
scribed in terms of a unitary transformation connecting the two input fields and 
the two output fields [^]. With regard to conditional measurement, it is conve- 
nient to regard the combined action on the signal of the beam splitter and the 
measuring instrument as the action of an optical two-port device. Following this 
concept, a non- unitary transformation operator in the Hilbert space of the signal 
field can be introduced which is independent of the signal quantum state. In 
this paper we present this operator in a closed form for arbitrary input quan- 
tum states of the reference mode and arbitrary quantum states measured in the 
output channel of that mode. 

The developed formalism generalizes and unifies previous approaches to the 
problem of conditional measurement on a beam splitter ^, |^, Q and enables us 
to calculate the conditional output states in a very straightforward way. To illus- 
trate the formalism, we consider the generation of Schrodinger-cat-like states from 
coherent and Fock states and give a brief extension to the generation of multiple 
Schrodinger-cat-like states. Originally introduced for probing the foundations of 
quantum mechanics, quantum superposition states of Schrodinger-cat-type (for 
a review, see [0) have recently been suggested to be applied as logical qubits in 
quantum computing P . 



The paper is organized as follows. Section |^ introduces the basic-theoretical 
concept and presents the non-unitary transformation operator. In Section ^ the 
formalism is applied to the generation of Schrodinger-cat-like states, and a sum- 
mary is given in Section 3. 
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2 The conditional beam splitter operator 



Let us consider an experimental setup of the type shown in figure |T]. A signal 
mode (index 1) is mixed with a reference mode (index 2) at a beam splitter, and 
a measurement (device M) is performed on the output reference mode. Since 
the two output modes are entangled in general, the measurement influences the 
output signal mode as well. The action of a beam splitter can be described by a 
unitary operator U connecting the input and output states according to 



where IT 
with 



|^out) = f>|^in), (1) 

fj — gl(</'T+</3jj)L3g2il9L2gl('/'T-'/5jj)i3 ^2) 

1 - 1 

L2 = Y^{a[a2 - a^ai), L3 = -(a|ai - 4«2) , (3) 

and the complex transmittance T and reflectance R of the beam splitter are 
deflned by 

T = e''^^ cos R = e''^« sin (4) 

Now let us assume that is the positive operator valued measure (POVM) 
that is realized by the measuring device M, with 

fi(/)>o, En(0 = i (5) 



(for POVM, see, e.g., [15, 16|). When the measurement on the output reference 
mode yields the result /, then the reduced state of the output signal mode becomes 



^?outi 



^outn(0 



where 



&ut = f>|^in)(^in|^7^ (6) 



pil) = (n(Z)) = TriTr2[foutn(/)J (7) 

is the probability of obtaining the result I. In particular, when n(/) projects onto 
a pure state \l) = |\E'out2)! i-e., 

n(0 = l^out2)(^out2l, (8) 

and the (pure) input state can be decomposed as 

|^i„) = |^in,)|^in2), (9) 

then combination of equations (^ - (^ yields 

&Utl = |^outi)(^outi|, (10) 
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with 



where 



I^l^in.) 



Y = (^outjt/l^ina), 



(11) 



(12) 



is the non-unitary conditional beam sphtter operator Y defined in the signal-mode 
Hilbert space, the expectation value of Y'^Y being the probability of obtaining 
the state I^E'outi), 



M^outJ = \\Y\^irn)f = i^inAY^Yl^irn)- 



(13) 



In order to determine the non-unitary transformation operator Y, let us first 
consider reference modes that are prepared in displaced Fock states (for displaced 
Fock states, see, e.g., and references therein). 



l^ina) = I^2(a)|m)2, l^out^) = ^2(/?)|^)5 



(14) 



with -D(a) =exp(aa^ — a*d) being the coherent displacement operator. After a 
lengthy but straightforward calculation (see [Appendix 7\|) we find that 



Y = Di\ — I 2{n\U\m}2Di\ 



R* 



R* 



(15) 



where 



{n\U\m)- 



R'^i-R* 



m\n\ 



J s 



(16) 



Here, the notation {■ ■ -js introduces s-ordering (for s-ordering, see 1T8|), with 
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\R\ 



1. 



(17) 



Note that s > 1 . Applying equation ( |A12| ) for t = 1 and using the formulas given 
in Appendix A in Ref. the s-ordered operator product in equation (|I^) can 
be rewritten as 



s + l 
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p^n—m p(n— m,n— 7i) 
(ci^)™~" p(m— ra,n— n) 



s-3 



s+l 
s-3 



s + l 



if m < n. 



if m > n. 



where P^'^\z) is the Jacobi polynomial. It should be mentioned that when a = 
(3 = 0, then the operator Y in equation ( plSj) realizes the transformation to the 
Jacobi- Polynomial states in Ref. H. 
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Now, the generalization of the formahsm to arbitrary (pure) quantum states 
of the reference modes, 

l^in^) = ^(4)|0)2, l^out^) = ^(4)10)2, (19) 

is straightforward. From equations (|T^ and (^) and apphcation of the relation 
(|16D Y is obtained to be 

R 



Y = !^F{Rd\)&^--d\^]T'^\ (20) 



The result reveals, that (up to the operator T"^) the non-unitary conditional 
beam splitter operator Y is nothing but the operator product FG^ in s order, 
with s from equation ([T7|). Although equation ( pO] ) already covers the general 



case, it may be useful, for practical reasons, to consider explicitly coherently 
displaced quantum states of the reference modes, that is 

l^in^) = ^2(a)F(a^)|0)2, l^out^) =^2(/?)G(a^)|0)2. (21) 
In close analogy to the derivation of equations (plSj) and (^) we find that 

>-«.(^){.(-l)c^'(-^.0}/--'(^)^ <-) 

which shows that displacing the quantum states of the reference modes always 
leads to displaced quantum states of the signal modes. 

From equation ( P^D it is easily seen that the ordering procedure can be omitted 
if at least one of the reference modes is prepared in a coherent state (the vacuum 
included), i.e., F = I oi G = I. This is the case when, e.g., |\E'm2) = |0)2 and 
|^out2) = 1^)2 (or |\E'm2) = 1^)2 and |\E'out2) = |0)2), which leads to preparation of 
the output signal mode in a photon-subtracted (or photon-added) state 0, |^. 
Note that preparation of the output reference mode in a coherent state can be 
realized in eight-port homodyne detection or heterodyne detection. Further it 
should be mentioned that I^E^im) and |^E'in2) can be interchanged if T and R are 
replaced with iR and iT, respectively, because of the symmetry of the beam 
splitter transformation \T^. The operator Y then transforms (up to a global phase 



factor) the state \'^m2) into the state l^&outi)- In this way the s-ordering procedure 
can also be circumvented when the signal mode is prepared in a coherent state. 

In general, the POVM realized by the measurement apparatus does not project 
onto a pure state, and equation (||) must be replaced with 

n(/) = Yl P(/|^OUt2)I^OUt2)(*OUt2l , (23) 
*out2 

where p(Z|\E'out2) is the probability of obtaining the measurement outcome / under 
the condition that the output reference mode is prepared in the state |\E'out2)- A 
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typical example is direct photon counting with quantum efficiency t] less than 
unity, 



(24) 



Further, the input reference mode may also be prepared in a mixed quantum 

st ate, 

ftn2 = EP(^in2)l^in2)(*inJ, (25) 

SO that for an input quantum state 

^in = ftm ® ftna (26) 

the quantum state of the output signal mode now reads 

Qout, = 4y E P(*in2) E p(/|^outjr^i,,Ft (27) 
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in place of (|TT]), where Y is defined by equation (|T2|). In equation ( p7D 



p(/) = E P(*in2) E p(/|*out2)Tri(rftn,V't) . (28) 

*in2 *out2 

is the probability of obtaining the measurement outcome I, i.e., the probability 
of producing the output signal-mode quantum state ^outi- 



3 Creation of Schrodinger- cat- like states 

In order to illustrate the theory, let us consider the generation of Schrodinger- 
cat-like states. A possible way is to prepare the input signal mode in a squeezed 
vacuum, combine it with an input reference mode prepared in a Fock state (in- 
cluding the vacuum state) and perform a photon-number measurement on the 
output reference mode 0, |^. Here we present two alternative schemes that are 
only based on Fock states and coherent states. Let us ffist consider a scheme 
[figure @(a)] that uses a Fock state source and displaced photon-number measure- 
ment, 

l*in,) = \n), (29) 



\%n,) = \0), l^out,) = Dmn). (30) 

For notational convenience we omit the subscripts 1 and 2 introduced in Section 
1^ in order to distinguish between the channels. Application of equations ([TT| ) and 



5 



J5[) then yields (|r|2 = 0.5) 



l^out,) =: Ix::-^) = -^(a-/5r(at + /5T|0) 

nW M 



fc=0 



N-^'^h^[(3D\(3)a^D{(3)\ |0) (31) 



[/3 = /3'e'('^'^+'^^+'^^; L„(z), Laguerre polynomial; L5^(2;), associated Laguerre poly- 
nomial], where 

n \fj\2k 



2-n,-W N = E ^ Lt,(|/?|f = ^n, /?'), (32) 

fc=0 

p(n, /?') being the probability of generating the state |Xn'^) given /3' [see equa- 
tion (pr3|)]. Equation (^) reveals that the quantum state of output signal mode is 
obtained by applying an operator Laguerre polynomial L„ on the vacuum state.Q 
As it is seen from the in figure ^ plotted contours of the Husimi function 

Qia) = -|(a|x::'^)r = 4F|Ln[/5(«*+/31l'e-l"^ (33) 
n TTiV 

the state \Xn'^) exhibits for = n/2 two well separated peaks in the phase 
space, which are approximately located at ±i/5 [figure 0(b)]. It should be noted 
that the distance of the peaks increases with the square root of the detected 
photons, which is analogous to the behaviour of the Schrodinger-cat-like states 
considered in [Q. The Husimi Q function, which can be directly measured, e.g., 
in eight-port homodyne detection or heterodyne detection, is a rather smeared 
phase-space function, because of the included vacuum noise. More details of the 
structure of the state can be inferred from the Wigner function 

1 r+oo 

W{x,p) = - / dye''^Py{x-y,0\xl'^){x':,^^\x + y,0), (34) 



TT J-oo 

where \x, 0) is the quadrature-component state 



\x,y,) = (7r)-3exp(-ix2) E ^=H,(x)|A;) (35) 

fc=o v2'=fc! 

for y9 = [Hfc(2;), Hermite polynomial]. Using equations (|31| ) and (|35|) and applying 
standard formulas |TD|, the integral in equation (|3^) can be calculated to obtain 



+ E ^i-km')K-m{wn~'^{\z?)fi\ - PT^-^, — (36) 



^Note that the states Ix"''') are different from the Laguerre polynomial states introduced in 
Ref. pOl]. The latter can be produced in the scheme studied in Ref. 0. 



[z = a/2(x + ip)]- The Wigner function is a measurable quantity which however 
does not correspond to a POVM. Let us therefore also consider the quadrature- 
component probability distributions 



\{x,V\Xn'^)\' 
k=0 



k\ 



(37) 



which also contain all available information on the quantum state and can be 
directly measured in four-port homodyne detection. Plots of the Wigner function 



and the quadrature-component distributions of the state \Xn ) •\/^/2 are 

presented in figures § and ^ respectively. They clearly reveal the typical features 
of Schrodinger-cat-like states. The probability of producing the state for various 
values of n is shown in figure ^. 

Let us compare the scheme in figure |^(a) with the somewhat modified scheme 
in figure 0(b). In the latter scheme the input signal mode that is prepared in a 
coherent state is mixed with an input reference mode that is prepared in a Fock 
state, and a photon-number measurement is performed on the output reference 
mode, 



(3^ 



1^ 



out 2 



) = 1^)- 



(39) 



The scheme realizes the generation of a special Jacobi-polynomial coherent state 



l^outi) ~ 



{(a+)"a"}s=3r"|/5/r) ~ P^"-")(0) 



(40) 



(|Tp = 0.5). Here we have used equations (0), (0) and (^) and the relation 
(fc) ~ {a^Ya^ /k\. Comparing equations ( PT| ) and (^), we see that the states pro- 
duced in the schemes in figures ^(a) and ^(b) differ in a coherent displacement. 
Obviously, the coherently displaced state shows again for = n/2 the typical 
features of a Schrodinger-cat-like state, the phase difference between the compo- 
nent states being 7r/2. The probability of producing the state is the same as in 
equation (^). 

It is worth noting that replacing in equation (^TJ) [a — (3) with [a — (3)'^ yields 
the state 



(at-/3f (at+/3f |0) = -^[(at)^-(r) 



|0), (41) 
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the properties of which are similar to those of the state |Xn^) equation (pi]). 
States of the type (0) may be produced, e.g., by two displaced n-photon- 
additions (for such schemes, see 111]). The generalization to k displaced n-photon- 
additions is straightforward. 



|0), 



States of this type can be regarded, for appropriately chosen /?, as multiple 
Schrodinger-cat-like states |^|, as it can be seen from figure |^, in which the 
Husimi function 



2fc(n-i) 



(42) 
(43) 



Q{a) 



1 



2n 



■lap 



(44) 



of such a state is plotted. 



4 Conclusions 

In this paper we have studied conditional quantum-state generation at a beam 
splitter, regarding the apparatus as an effective two-port device, which is de- 
scribed by a non-unitary beam splitter operator Y that acts in the Hilbert space 
of the signal mode only. We have presented the Y operator for arbitrary quantum 
states \'^m2) and |\I/out2) of the reference mode as an s-ordered operator product 
of those operators which generate |\I^in2) and (\l'out2l from the vacuum, s being 
entirely determined by the absolute value of the reflectance (or transmittance) of 
the beam splitter. We have further given a generalization to the case of mixed 
input states and non-perfect measurement. 

The formalism is a generalization of previously obtained results, which are 
easily recovered by appropriate specification of the non-unitary transformation 
operator. Given the quantum state of the input signal, the conditional quantum 
state of the output signal can be calculated in a very straightforward way for 
arbitrary input reference states and arbitrary measurement-assistant projection. 
The probability of producing the conditional quantum state of the output signal 
is simply given by its norm. In order to illustrate the theory, we have considered 
two schemes for generating Schrodinger-cat-like states and addressed the problem 
of generating multiple Schrodinger-cat-like states. Finally, it should be pointed 
out that the formalism developed is also suitable for studying state manipulation 
and state engineering via conditional measurement at multiports that are built 
up by beam splitters, such as chains of beam splitters. 
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Appendix A Derivation of Equation 

From equations (|12D and (0) the operator Y is defined by 

Y= ,{n\Dl{f3)UD2{a)\m)2, (Al) 
where U, equation (H), can be rewritten as 



(A2) 



We insert equation (^^) into equation (|A1|) and calculate the channel-2 matrix 
element as 



where we have used the relation 

and |n) = l/v^(a^)"|0). Applying the relations 

after straightforward algebra we derive 



y= ^ EE 



-a +-ar 



n—l 



2{m 



Using the ordering formula [18 



min{m,n} 

E 

1=0 



I ) {n-l)\ 

together with equation and the relation 

T"|a) = e-l-l'(i-l^l')/2|Ta) 

we have 



,a/3VTg-(|al2 + l/3|2)/2^^ ^ 



^ n min{fc,Z} /^^^ ^^^^ 



A; / U 



^m\n\ fc=o z=o j=o 



R 



m—k 



*\k~j 



(A3) 



(A4) 



(A5) 



^(a5)^T-"2|a)2e(^/^)""l (A6) 



(A7) 



(A8) 



Changing the summation indices, some of the finite summations may be per- 
formed to obtain 

rphi-m min{m,7i} / \ / \ 

Y = /Tg-(|a|2 + |/3|2)/2 ( ) ( ^ 1 fcl 

VM. to \kj\kj ■ 

X e-^*''*'^ (^-R*a, + ^-^y (i?aI-Ta*+/3*)"^"'e(^/^)"'^l (AlO) 

Applying again the relations ( |A^ ) and ( |A5D , after some calculation we obtain 



-l-(2|i?|-2-l)''= 



T J ^ ' ^ ^ ^' \ '\ R* 

Making use of the standard ordering formula 



T""^ ^ /a-r/?\ f™"^""'"^ (m\(n^ 
yi{-\R\''f (-—\ {RT)'^-^a'l~\a\)'^-^T^AbA ^~^*'^ \ . (All) 



min{m,n} / \ / \ ,x \ k 

' m\ n\ f t — 



k=0 



for t = — 1, we eventually arrive at equation (^) (for more details, see ||22|). 
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Figure 1: After mixing a signal mode prepared in the state l^&mi) with a reference 
mode prepared in the state |\I/in2)) the output signal mode collapses to the state 
l^outi), if the measuring instrument M projects the output reference mode onto 
the state |\l/out2)- 
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Figure 2: (a) A Fock state \n) is fed into a balanced beam splitter. When the 
measuring instrument M projects the output reference mode onto a displaced 
Fock state D{(3')\n), then for = = n/2 the output signal mode collapses 
to the Schrodinger-cat-like state \Xn^)i equation (pT|). 

(b) A coherent state \I3/T) and a Fock state \n) are fed into a balanced beam 
splitter. When the measuring instrument M projects the output reference mode 
onto the same Fock state |n), then the output signal mode collapses to the state 
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Figure 3: Contour plots of the Husimi function Q{a), equation (^), of the state 
Ix"''^), equation (|31|), for n = 10 and (a) /3 = yn/lO, (b) l3=^n/2, and (c) (3 = 
VS. 
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Figure 4: The Wigner function W{x,p), equat 
equation (^), is shown for n = 10 and /? = \Jn/2. 
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Figure 5: The quadrature-component probability distributions p{x, cp), equation 
(p?!), of the state |Xn'^); equation pl|), are shown for n = 10 and /? = Jn/2. 
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Figure 6: The probability p{n,i3'), equation (|32D , of producing the state \Xn^)y 
equation (pl|), is shown as a function of n for = = n/2. 
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